How to count nucleon pairs? 
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Abstract 



Within the context of the isovector-pairing S0(5) model, three methods measur- 
ing numbers of different kinds of nucleon pairs are discussed. Though methods do 
\ not give the same results, the odd-even staggering in pair numbers in even-even and 

odd-odd N = Z nuclei and the reduction of the np-pair number with increasing 
is observed in all three procedures. 
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^ Structures comprised of nucleon pairs are influential in describing nuclear 

states. The most important pairing degree of freedom is the monopole pair 
with L = orbital angular momentum. In medium-heavy and heavy nuclei, 
^ I the monopole pairing pertains to pairs of like nucleons as valence protons and 

Ch ■ neutrons occupy different shells and cannot couple to L = 0. Today, consid- 

erable experimental activity is aimed at study of proton-rich nuclei. In those, 
^ ! the neutron-proton-pairing mode becomes notable. 

H ; 

In a recent paper fl]], a method has been suggested to quantify effects of differ- 
ent pairing modes and measure pair numbers in the nuclear-wave function. We 
should emphasize that the pair numbers are not directly observable quanti- 
ties. Even in simple models of pairing, there is no operator that counts nucleon 
pairs. If an exact solution of the problem is available, one need not to operate 
with the pair numbers. Their quantification can, however, make the problem 
more transparent. Also in complex situation, the notion of the significance of 
different pairs may be helpful in constructing approximate solutions and/or 
devising model descriptions. 

In the present study, we examine the measure of pair numbers of Ref. ^ 
and suggest alternative procedures. As in Ref. [|I|, we discuss the problem 
within the context of a simple model of monopole-isovector pairing with the 
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underlying 0(5) algebraic structure p[. The model reflects pairing features of 
realistic shell-model calculations 0. 



In the monopole-isovector-pairing model, the set of pair creation operators 

Sl = ^i:[a'^-a^''-]l , (1) 

the conjugate pair annihilation operators 5*,^, the isospin operator T, and the 
total number operator N closes the 0(5) algebra. The degeneracy of the shell 
is 4f2 = 2^{2j + 1). The problem of the S0(5) pairing Hamiltonian 

H = -GY.SlS,, (2) 

V 

can be solved algebraically. We shall consider the seniority zero eigenstates 
\MTTz) of H labelled by half the number of nucleons A/", isospin T, and isospin 
z-projection T^. 

The basic building blocks of the SO (5) model space are the nn-, pp-, and np- 
pair creation operators. We shall examine possible procedures which would 
estimate the average numbers (A/'nn), (A/'pp), and (A^np) of the respective pairs 
in a given state. 

Within the 0(5) algebra, there are no operators which could be considered as 
the number operators of particular pairs. There are, however, two operators 
which relate the different pair numbers. Namely, the sum of the pair numbers 
should give half the number of nucleons 

(ATnn) + (ATpp) + (Ar„p) = M (3) 

and the difference between nn- and pp-pair numbers should link with the 
charge of the state 

{Unn) - (ATpp) = . (4) 

In Ref. the operator 

A/'nn = yiSl (5) 

and the analogously defined operators N^^ and AT^^ have been suggested as 
measures of the numbers of different pairs. We denote this method by I. In 
this method, the pair numbers are associated with the contributions of the 
respective pairing terms in the Hamiltonian (^ to the total energy and the 
relations can be obtained 

{K.) + «) + {K) = m - ^) - ^T(T + 1) (6) 
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«)-«)=T.(l-:^) (7) 



3(2T-l)(2T + 3)^^ 1] j(2A/+3j 

.l(Ar + T + l)(Ar-T)] . 



For the large shell degeneracy Q or more strictly speaking for -C O, the pair 
operators exhibit boson-like behaviour and the pair numbers in method I obey 
Eqs. @ and (Q). For nuclei around the middle of the shell, however, serious 
discrepancy between and and the natural constraints (H) and on 
pair numbers may occur. One can remedy the total pair number condition by 
scaling simply pair numbers in method I (and we shall use this procedure in 
the following). Anyhow, the charge constraint (|) will not be fulfilled. 

This deficiency puts some question on the method I as a measure of differ- 
ent pair numbers. Certainly, operators quantify contributions of different 
pairing terms in Hamiltonian (|]). In many-pair states, however, those are in- 
fluenced by the Pauli principle and the relation to the pair numbers is not 
immediate and even may be misleading. 

We shall propose an alternative method II for obtaining the pair numbers that 
employs directly the pair basis. In that the states are constructed as 

lAT^^ATppAT^p) = aSf^-Sif^Si^-^lO) , (9) 

where a normalizes the state (P). The use of the pair basis in calculating the 
average pair numbers is not straightforward. Namely, the pair basis is not 
orthogonal and for A/" > is even overcomplete. Nevertheless, the squared 
overlap of the pair state @ with a wave function lip) gives certainly a prob- 
ability of the occurrence of the state with given pair numbers. We are thus 
tempted to define the average pair numbers as 

(Ki) = Cj^KMnMpP^nMl' (10) 

with the scaling factor 

C = El(-^n„ArppAr„p|^)|2 

and similarly for (A/'pp) and (A/'jJp). Note that C scales the definition ( |I0| ) so 
that Eqs. (|^) and are obeyed. In the limit of large Q, pair numbers from 
method II agree with those from method I. 
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The complications in the preceding two definitions (j^) and ([lOD of pair num- 
bers are connected with the fact that we intend to implement the boson-like 
behaviour to pair operators that do not obey the boson commutation relations 
exactly unless Q is large. This mirrors in the insufficiency of the method I to 
fulfil the charge constraint (^) and in the necessity to deal with the nonorthog- 
onal basis in the method II. 

There exists a possibility to rephrase the fermion problem in the boson lan- 
guage. H For the bifermion 0(5) algebra, the Dyson boson realization has been 
constructed in Ref. In our specific isovector-pairing case, the boson space 
is formed from the scalar-isovector boson s^,. The fermion-pair operators are 
mapped on the boson operators as 

S,^s,. (11) 

Using the above relations, the boson image of the Hamiltonian H is con- 
structed. When Hb is diagonalized in the ideal (Fock) boson space, all the 
eigenvalues and boson images of eigenstates of the original fermion problem 
are provided. The dimension of the ideal boson space can be larger that that 
of the fermion task. For the SO (5) Hamiltonian, spurious states occur in the 
boson picture for J\f > Q. They, however, clearly separate from the physical 
solutions for example by the isospin value, i 

We may make use of the boson formulation of the fermion problem to get 
the boson-like fermion-pair numbers. The bosons are images of the fermion 
pairs. In the boson space, the boson number operators n are well defined. We 
identify the average boson numbers with another measure (method III) of the 
pair numbers 

(Kn) = (K) (12) 

and similarly for {Af™)a.iad (A/'^p). Here, the boson eigenvectors are denoted 
by the round brackets. □ Such a definition of pair numbers obeys Eqs. (0) and 
(11). Explicitly, one gets 



^ For a detailed review on boson realizations of fermion Lie algebras see . 

^ The boson mapping may provide a simple technique to solve the fermion problem. 
In fact, we use it to obtain the pair numbers in both methods I and II . 

^ Note that the Dyson boson mapping is a nonunitary mapping so that the boson 
Hamiltonian may generally be nonhermitian ( this is not the case for the Hamilto- 
nian ^) with different left and right boson eigenvectors. 
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Fig. 1. The pair numbers for the S0(5) model ground state in = Z nuclei. The 
filled circles, triangles, and hollow circles refer to methods I, II, and III, respectively. 



Alternatively, the boson definition (12) of pair numbers is obtained from the 
method I employing the pairing force when the 1/VL corrections are neglected, i 
One cannot, however, deduce that the definition (|T2|) is inferior to that of (^. 
The 1/Vt corrections reflect the Pauli principle and deviations of the pair be- 
haviour from the boson one. On the other hand, we would like to understand 
the pair numbers which are the boson-like quantities. Simple relations of Eqs. 
d^) and @ for the pair numbers should hold which are in the method I dete- 
riorated just by the 1/VL terms. 

As it is seen from Fig. |l], all three methods give very similar quantities for 
the average pair numbers in the ground states of N = Z fp shell nuclei. The 
odd-even staggering observed in Ref. [0 is thus confirmed also by another two 
measures of the numbers of different pairs, i 

In Fig. D on the left, the pair numbers for even-even Fe isotopes are shown. 
Here, the three methods differ more pronouncedly than for the case of Fig. 



^In the present S0(5) case, we know the results in method I analytically and can 
perform the Q ^ oo limit without resorting to the boson formulation. In the general 
case, the diagonalization in the boson space is necessary. 

^ Note that the proximity of the results of the methods I and III does not mean 
that the l/Cl corrections in are small. We remind that the pair numbers from 
the definition dsl) are scaled to sum to the total number of pairs. 
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Fig. 2. The pair numbers for the S0(5) model ground state in even-even Fe and 
odd-odd Co nuclei. The filled circles, triangles, and hollow circles refer to methods 
I, II, and III, respectively. The full lines connect (AAnn) points, the dash-dotted lines 
connect {^fpp) points, and the dashed lines connect {M^^^ points. 

|1]. Nevertheless, the alternative methods confirm a decrease of the number of 
np pairs with increasing noticed in Ref. f^j- In the method I, numbers of 
nn and pp pairs increase with increasing and significantly break down the 
charge condition (1). On the other hand, in methods II and III that fulfil the 
charge constraint, a rough constancy of the number of pp pairs and linear 
increase of the number of nn pairs is observed when the neutrons are added to 
the nucleus. In Fig. ^ on the right, the odd-odd Co isotopes exhibit a similar 
behaviour with increased importance of the np pair as compared to even-even 
case. 

When Tz is increasing, the decrease of the number of np pairs suggests a 
diminishing role of the np pairing. It would be important to be able to judge 
at what value of (A^ — Z) one can neglect the np-pairing degree of freedom 
and approximate the ground state as a product of nn and pp paired states. 
The exactly solvable SO (5) model makes possible such a discussion. The exact 
ground-state wave function of the Hamiltonian for the even-even nucleus 
is written up to the normalization factor as 

whereas for the odd-odd nucleus we have 

\AfT = T, + lT,) ~ {S^ ■ S^)(^-^^-i)/'5pS(l|0) . (15) 
The important building block of the ground-state wave function is the two-pair 
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Fig. 3. In the upper part, the exact ground-state energy (in units of G) (hohow 
circles) and the mean value of H in the approximate state with no np correlations 
(filled circles) are given. In the lower part, the squared overlap of the exact and 
approximate wave functions is shown. 

T = structure 

-S^ = Sd^d - 2Slsl^ . (16) 

From that the T = even-even core is constructed to which an appropriate 
number of nn pairs and one np pair for the odd-odd case is added. 

We shall approximate the SO (5) ground state by discarding the np part in 
the above two-pair T = structure (0). The isospin is then mixed in the 
approximate solution. For the even-even nucleus, we have 

\XT, , no np correlations) ~ 5t(Ar+T.)/2^KAr-T.)/2|Q^ _ ^^^^ 
Of course, at least one np pair must be present for the odd-odd nucleus 

|AAT, , no np correlations) ~ Sj^^+'^^-'^/'slf ""^^"'^/'sdlO) . (18) 

In the upper part of Fig.^, the exact ground-state energy is compared with 
the mean value of H in the approximate state with no np correlations. In 
the lower part of Fig.^, the squared overlap of the exact and approximate 
wave functions is shown. Fig.Q corresponds to the findings inferred from the 
analysis of pair numbers. The diminishing role of the np pair with increasing 
Tz is particularly exhibited in the even-even case. In the odd-odd nucleus, the 
mandatory np pair stabilizes the np-pairing degree of freedom and there is 
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a little overlap between the exact state and the wave function with the np 
correlations neglected for small values of T^. Here, however, the increase of 
the overlap is more abrupt than in the even-even case. For (A^ — Z) = 10 both 
in the even-even and odd-odd case, the exact wave function overlaps with the 
wave function with no np correlations by more than 90%. 

Note also that the exact and approximate ground-state energies might not 
differ much even if the overlap is small. The leading term in the ground-state 
energy is given as —GQf/ irrespective of the actual form of the wave function. 
The binding energies alone are not thus best quantities to judge about the 
np correlations and about the quality of an approximate wave function. One 
should resort to difference filters, such as that of Ref. in which the leading 
term contribution is eliminated. 



If the isospin symmetry is violated in the SO (5) Hamiltonian and the np- 
pairing force is set to zero, one would expect that the number of np pairs 
should go to zero. However, the methods I and II give a nonzero average 
number of np pairs in this case. Similarly, in the 0(8) space including the 
isoscalar pair 0], the number of isoscalar pairs is nonzero for the Hamiltonian 
with the isoscalar pairing switched off in methods analogous to and (TU). 
On the other hand, the boson method HI of counting pairs gives a null number 
for the pairing mode whose pairing term is not present in the Hamiltonian. 



To conclude, we stress again that counting the number of pairs in a fermion 
state is not well and completely determined task. There is no fermion operator 
which could be identified with the number operator of a particular pair. Nev- 
ertheless, three possibilities have been investigated in the present paper how 
to get quantities that may be related to the pair numbers. Though methods 
give different results, two important conclusions of Ref. remain unchanged. 
Namely, the odd-even staggering in pair numbers in even-even and odd-odd 
N = Z nuclei and the reduction of the np-pair number with increasing are 
observed in all three procedures. The boson counting method seems to have 
certain advantages as it obeys the natural relations (|^) and for pair num- 
bers and gives a zero pair number for the particular pairing mode switched off 
in the Hamiltonian. 
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